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The  overall  objective  of  the  work  done  under  this  contract 
is  to  determine  the  feasibility  of  using  solitons,  rather  than 
sinusoidal  wavetrains,  as  the  basic  building  blocks  in  a theory  of 
water  waves.  The  goal  of  devising  such  a revised  description  is 
to  model  phenomena  that  occur  over  longer  time  scales  than  can  be 
adequately  described  in  terms  of  sinusoidal  wavetrains.  This 
report  summarizes  the  work  done  in  this  direction  during  the  first 
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The  work  done  under  this  contract  resulted  in  two  scientific 
papers,  both  of  which  were  co-authored  with  Mark  Ablowitz: 

(i)  "Long  internal  waves  in  fluids  of  great  depth," 
submitted  to  The  Physics  of  Fluids;  and 

(ii)  "On  the  evolution  of  packets  of  water  waves,"  sub- 
mitted to  the  Journal  of  Fluid  Mechanics. 

These  papers  are  include  in  this  report  an  Appendices  A and  B.. 

The  first  paper  discusses  three-dimensional  effects  on  the 
kind  of  internal  waves  first  discovered  by  Benjamin  (1967)  and 
Davis  and  Acrivos  (1967).  In  this  paper,  we  derive  the  ap- 
propriate higher-dimensional  generalization  of  Benjamin's  equa- 
tion describing  the  evolution  of  these  waves,  and  show  that  the 
(plane-wave)  solitons  are  not  unstable  to  long  transverse  per- 
turbations . 

The  second  paper  represents  the  major  effort  under  this 
contract.  In  it,  we  discuss  three-dimensional  effects  on  packets 
of  surface  waves,  including  both  long  waves  (governed  by  the 
Kor teweg-deVries  equation  in  two  dimensions)  and  short  waves 
(governed  by  the  nonlinear  Schrodinger  equation  in  two  dimensions) . 
In  each  case,  we  derive  the  appropriate  higher-dimensional  equa- 
tion, along  with  the  appropriate  boundary  conditions.  Then  using 
these  equations  we 

(i)  analyze  the  stability  of  (plane-wave)  solitons  to  long 
transverse  perturbations; 

(ii)  show  the-existence  of  "lumps,"  which  are  higher  dimen- 
sional analogues  of  solitons; 

(iii)  derive  conditions  under  which  a packet  of  capillary- 
type  waves  must  "focus"  at  a point  in  a finite  time 
(a  strong  nonlinear  instability); 

(iv)  investigate  the  suitability  of  Inverse  Scattering  Trans- 
forms for  these  higher  dimensional  equations;  and 

(v)  discuss  some  special  solutions  that  may  have  physical 
interest . 
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It  is  difficult  to  describe  the  outcome  of  these  analyses 
simply,  because  different  results  are  obtained  in  different 
ranges  of  the  dimensionless  fluid  depth  and  dimensionless  sur- 
face tension.  Certainly  the  evolution  of  wave  packets  in  three 
dimensions  involves  phenomena  that  cannot  occur  in  two  dimen- 
sions. More  detail  has  been  given  in  the  previous  quarterly 
progress  reports,  and  of  course  in  the  paper  itself  (Appendix  B) . 

The  work  in  these  two  papers  shows  the  direction  of  research 
that  seemed  most  fruitful  as  it  developed.  The  final  result, 
however,  is  not  obviously  related  to  the  set  of  problems  out- 
lined in  the  original  proposal  for  this  contract.  Therefore, 
it  may  be  worthwhile  to  comment  hereon  the  current  status  of 
the  problems  that  were  originally  outlined  almost  two  years 
ago. 

1.  It  was  proposed  to  analyze  the  resonant  interaction  of 
a long  internal  wave  with  a packet  of  short  surface  waves, 
originally  discussed  by  Phillips  (e.g.,  1974).  The  proposed 
analysis  would  employ  multiple  time-scales,  in  order  to  see  if 
the  evolution  could  be  described  in  terms  of  solitons.  The 
answer  is  affirmative,  but  no  work  on  that  problem  was  done  under 
this  contract,  primarily  because  the  analysis  had  already  been 
initiated  by  Larry  Redekopp  (private  communication) . His  work 
should  appear  shortly  as  a TRW  Report;  related  work  was  also  done 
by  Grimshaw  (1977)  . 

2.  It  was  proposed  to  analyze  the  oblique  interaction  of 
envelope  solitons,  to  see  whether  an  interaction  like  Miles' 

(1977)  three-wave  resonance  of  KdV  solitons  could  occur.  The 
answer  is  affirmative,  as  shown  by  Newell  and  Redekopp  (1977). 
However,  the  results  of  our  stability  analyses  are  relevant  here. 
The  fact  that  envelope  solitons  are  unstable  to  long  transverse 
perturbations,  whereas  KdV- type  solitons  ordinarily  are  not,  sug- 
gests that  the  oblique  interaction  of  envelope  solitons  has 
limited  physical  significance. 
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3.  It  was  proposed  to  reinterpret  field  measurements  of 
ocean  wave  spectra  in  terms  of  solitons.  But  these  spectra  are 
decidedly  three-dimensional  in  character.  The  differences  al- 
ready discovered  between  wave  evolution  in  two  and  three  dimen- 
sions suggest  that  (two-dimensional)  solitons  alone  will  not 
explain  these  spectra,  and  points  Lu  Lhe  need  for  more  work  on 
the  three-dimensional  problem. 

4.  It  was  proposed  to  study  the  stability  of  envelope  soli- 
tons to  transverse  perturbations.  This  analysis  was  completed 
under  this  contract.  Experiments  to  test  these  results  are  now 
being  set  up  by  Joe  Hammack. 

5.  It  was  proposed  to  study  the  viscous  decay  of  envelope 
solitons.  The  experimental  portion  of  this  problem  has  been 
completed;  development  of  an  appropriate  theory,  and  its  com- 
parison with  the  experimental  data  should  occur  within  the  next 
year . 

In  summary,  the  work  completed  this  year  on  the  evolution 
of  packets  of  water  waves  in  three  dimensions  indicates  that 
even  though  the  two-dimensional  theory  (involving  solitons) 
may  describe  well  the  essentially  two-dimensional  laboratory 
experiments  on  water  waves,  neither  is  necessarily  representative 
of  what  happens  in  the  three-dimensional  ocean.  We  have  begun 
to  understand  some  of  the  additional  phenomena  that  occur  in 
three  dimensions,  but  much  more  work  is  required. 
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APPENDIX  A 


LONG  INTERNAL  WAVES  IN  FLUIDS  OF  GREAT  DEPTH 


by 

M . 3 Ab 1 owi 1 2 ^ 
Program  in  Applied  Mathematics 
Princeton  University 
Princeton,  New  Jersey  08540 

and 


H.  Segur 

Aeronautical  Research  Associates 
of  Princeton,  Inc. 

50  Washington  Road 
P.O.  Box  2229 

Princeton,  New  Jersey  08540 


ABSTRACT 

An  equation  is  derived  that  governs  the  evolution  in 
two  spatial  dimensions  of  long  internal  waves  in  fluids 
of  great  depth.  The  equation  is  a natural  generalization 
of  Benjamin's  (1967)  one-dimensional  equation,  and  relates 
to  it  in  the  same  way  that  the  equation  of  Kadomtsev  and 
Petviashvili  relates  to  the  Korteweg-deVries  equation. 

The  stability  of  one -dimensional  solitons  with  respect  to 
long  transverse  disturbances  is  studied  in  the  context 
of  this  equation.  Solitons  are  found  to  be  unstable  with 
respect  to  such  perturbations  in  any  system  in  which  the 
phase  speed  is  a minimum  (rather  than  a maximum)  for  the 
longest  waves.  Internal  waves  do  not  have  this  property, 
and  are  not  unstable  with  respect  to  such  perturbations. 

♦Permanent  Address:  Department  of  Mathematics,  Clarkson 
College,  Potsdam,  New  York,  13676. 
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I.  INTRODUCTION 

It  is  known  that  the  evolution  of  long  waves  of 
moderate  amplitude  in  a nonlinear  system  depends  on  the 
form  of  the  linearized  dispersion  relation.  If,  in  the 
long  wave  limit  (k  -*■  0)  , the  linearized  dispersion  rela- 
tion is 

oj2  = c2  k2  - 2Bcok4  + 0(k6)  , (1.1) 

the  governing  equation  is  typically  the  Korteweg-deVries 
(KdV)  equation  [1,2,3]  (assuming  quadratic  nonlinearities), 


9u  , 9u  , 
r-  + au  r-  + 
9t  9x 


0 


If  the  dispersion  relation  has  the  form 


(1.2) 


u>2  = c2k2  - 2bcQk2|k|  + 0(k“)  , (1.3) 

then  (1.2)  is  replaced  by  the  equation  of  Benjamin  [A], 

+ au  |H+  b ^ H( u)  = 0 , (I-*) 

where 

H(£)  ’ ? / 3^f<y)dy 

is  the  Hilbert  operator.  [Note  the  sign  convention.]  Both 
equations  are  known  to  have  soliton  solutions  [5],  and 
these  agree  in  many  of  their  important  features: 
for  (1.2), 

u - (12K2  8/a)sech2 (K(x  - 4K2  Bt)}  (1.5) 
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and  for  (1.4),  with  cb  > 0 , 

u = ¥ {(I  (x  ' ct))2  + 1 • (1-6) 


Whereas  the  KdV  equation  can  be  solved  exactly  as  an 
initial  value  problem  [6,7],  a corresponding  method  of 
solution  for  (1.4)  has  not  yet  been  discovered.  There- 
fore, it  is  not  presently  known  how  different  the  general 
solutions  of  the  two  problems  really  are. 

If  the  long  waves  in  (1.2)  are  subjected  to  even 
longer  transverse  modulations,  Kadomtsev  and  Petviashvili 
[8]  reasoned  that  the  KdV  equation  should  be  generalized 
to : 


x (1.7) 


Moreover,  they  showed  that  the  soliton  solution,  (1.5), 
is  either  neutrally  stable  for  a limited  time  or  unstable, 
depending  on  whether  B is  positive  or  negative;  i.e., 
if  the  linearized  phase  speed  (w/k)  has  a maximum  at 
k = 0 (B  > 0)  , then  the  solitons  are  not  unstable.  For 
example,  this  is  the  situation  for  long  water  waves  without 
surface  tension. 

[For  B > 0 (the  stable  case) , no  information  travels 
faster  than  the  waves  described  by  (1.7),  and  it  is  ap- 
propriate to  take  xq  * + °°  . Further,  one  may  require 
that  u -►  0 as  x , but  then  one  cannot  require  that 


A 


-J 
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I 

| 
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u 


-*■  0 as  x = - 00  unless 


7*a2u 

J.  3y 2 


dx  = 0 


For 


6 < 0 , the  situation  is  reversed  and  one  should  take 

xq  = - 00  . Moreover,  it  is  worth  noting  that  when  the 

one-dimensional  solitons  are  unstable  (B  < 0)  , (1.7) 

admits  "lump'1  solutions,  which  are  localized  in  both  x 

'T  32u 

and  y and  do  have  J - — T dx  = 0 [9].  Their  interaction 

— OO  * 

behavior  suggests  that  lumps  are  two-dimensional  analogues 
of  solitons.] 

The  purpose  of  this  note  is  to  show  that,  with  regard 
to  very  long  transverse  modulations,  (1.2)  and  (1.4) 
behave  similarly.  Specifically,  we  will  show  that: 

(i)  the  appropriate  generalization  of  (1.4)  for  two-dimensional 


waves  is  (cq  > 0) 


9u  9u  , , 9 2 s o 

9t  aU  9x  + b 9^  H(u)  " T 


o 

/ 


9 2u 

w 


d£ 


(1.8) 


7 92u 

—5-  dx  = 0 ; and 
— 00  0 y 

(ii)  the  soliton  solution,  (1.6),  is  either  neutrally 
stable  for  a limited  time  or  unstable,  depending  on 
whether  b is  positive  or  negative.  Again,  if  the 
linearized  phase  speed  has  a maximum  at  k = 0 , then 
the  solitons  are  not  unstable.  This  is  the  situation  for 
long  internal  waves,  traveling  along  the  thermocline  in 
the  deep  ocean  [4,10,11]. 

In  Section  III,  we  show  how  to  derive  (1.8)  in  the  context 
of  internal  waves  on  an  ocean  thermocline.  This  derivation 
is  based  on  the  work  of  Ono  [11],  but  we  have  found  his  analy- 
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sis  to  be  incomplete.  In  particular,  the  matching  section  of 
his  analysis  applies  only  if  the  density  gradient  is  con- 
tinuous, although  no  such  restriction  is  stated  (or  even 
suggested)  in  [11].  In  order  to  complete  this  analysis, 
therefore,  we  show  how  to  match  correctly  both  where  the 
density  gradient  is  continuous,  and  where  either  the  density 
or  its  gradient  is  discontinuous. 

Finally,  we  note  that  if  one  is  simply  interested 
in  the  structure  of  the  underlying  equation,  such  a derivation 
is  unnecessary,  because  one  can  see  directly  from  the  linear 
dispersion  relation,  (1.3),  that  (1.8)  is  the  appropriate 
generalization  of  (1.4).  In  two  dimensions,  we  must  interpret 

k2  = k2  + k2 
x y 


in  (1.3).  In  order  that  the  effects  of  y-variations  enter 
at  the  same  order  as  the  linear  dispersion,  we  require 
that 


k’  = O(k^) 


<<  1 


(1.9) 


Then  (1.3)  becomes 


2 

U)  C 


■k’fl+(M-  2 — Ik  ! 

° xl  Vv  co  x J 


and  for  waves  in  one  direction,  we  obtain 


k u)  ~ c k2  - bk2  |k  | + k2 
x ox  x 1 x 1 2 y 


n .10) 


This  corresponds  to  a certain  linear  operator  acting  on 

exp  (i(k  x + k y - cot)}  . After  the  Galilean  transforma- 
x y 
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tion 

S = 

X - c t , 
o ’ 

t = t , this  operator  becomes 

<kx  + 

-i9x 

co  -*■  i3t 

► 

, sgn  (kx)  -iK) 

Lv  ’ ( 

9v,  , 8 2 u , A , Co  9 2 v n 

. d t + b h2H(v);  + t 37  = 0 

Integrating  once  gives  the  linearized  version  of  (1.8) 
and  the  analysis  which  yields  (1.4)  shows  that  the  non- 
linear term  also  enters  at  this  order. 
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II . STABILITY  OF  THE  SOLITOM 

In  order  to  test  the  stability  of  (1.6)  to  weak 
transverse  perturbations,  it  is  convenient  to  rescale 
y so  that  (1.8)  becomes 

x 

o 

dU  , &U  j.  u d2H(u)  x 2 I 9*U  j„i  _ A 

t—  + au  -r — ■+■  b — ~ — 7 — - - 4 I -5— r dx  = U , 

dt  9x  9x  I 9y 

x <5<<1  . (2.1) 

We  define  a slow  time  scale,  and  a spatial  scale  that 
travels  with  the  wave: 

T - fit 


z = x - f cdt'  - 6 (T , y , 6 ) 


(2.2) 


and  we  assume 


u = u (z,T,y,£)  + 6u^(z,T,y,6) 


(2.3) 


Then , 


9u  - 9u  , x 9u 
9t  = "C  91  + 6 9T 


9u  _ „ 9u  , 9u 

9y  r 9z  dy 


Substituting  into  (2.1)  yields 
9u  9u  ~2 

-jr^+au  j^+b4  H(U  ) 

9Z  09z  dZ  O 


c = c(T)  + 6 


p = - 9y  • 


[9u  9u,  a a2 

YT  - E YT  + a 91  (uoul)  + b 9^H(u1^| 

• fi2[?r  + aui  yt  + p2  yt  + 2p  9>r  + !y  uo 

f°  1 

J Ty£  dz’J  - °<63> 


0 (6 3 ) 


(2.4) 


(2.5) 
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To  leading  order,  (2.5)  is 


ctU  3u  ,2 

-c  + au  -r — + b H (u  ) = 0 , (2.6) 

3z  o 3z  dz  o 


and  the  solution  of  interest  is 


4c 

a 


u (z;  c)  TfzT"2- 
0 (¥)  + 1 


cb  > 0 


(2.7) 


At  higher  order,  the  equations  all  take  the  form 


Lf  = F , 


(2.3) 


where 


Lf  ’ M + * h <“of>  + b h*  H(f> 


The  usual  secularity  condition  is  found  by  multiplying 
(2.8)  by  uq(z)  « integrating  on  (-°°,°°)  , and  using  (2.6) 


J uoFdz  ■ [-r 


f - bH(u  ) 
o 


M 1 

3z 


Inserting  (2.7),  this  can  be  written  as 


00  <* 

/i7j  [4b2  d f z2f  \ ! 

uoFd2  ITT  3 z{  Z,  ‘ (|JZ)  J 


Thus,  the  solution  of  (2.8)  and  its  derivative  is  bounded 
for  all  z only  if 


oo 

/ 


u Fdz  * 0 
o 


(2.9) 


— OO 
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This  provides  our  basic  secularity  condition.  We  also 
note  that  (2.9)  will  not  yield  a uniformly  valid  solution 
on  | x | < * , but  does  give  the  correct  evolution  following 
the  soliton.  This  is  also  true  for  (1.7). 

At  0(6)  , (2.5)  yields 

3u 

Lui  = - (2.10) 

Applying  (2.9)  requires  that 

|f  - 0(6)  ; c = c + 6 |f  (T, y , 6)  (2.11) 


Hence,  the  first  nontrivial  result  comes  by  taking 
U1  = ^u2  ' The  9uesti-on  of  stability  is  determined  at 
0(62)  , where  (2.5)  becomes 

o 


Lu, 


3u 
c 

3T 


3uc 

3z 


- 2p 


3u  - 

- i£  u 
3y  3y  o 


J 


32u 

c 


dz 


(2.12) 


Applying  (2.9)  to  (2.12),  and  using  3c/3y  = 0(6)  , yields 
an  expression  which  reduces  to 

1 32G  . - 326  n 

■ITF+C  • (2-13) 

If  b < 0 , then  c < 0 (from  (2.7)),  and  (2.13)  shows 
that  arbitrarily  small  y-variations  in  the  'hase 
of  the  soliton  grow  rapidly;  i.e.,  the 

soliton  is  unstable  with  respect  to  such  transverse  per- 
turbations. In  fact,  (2.13)  is  elliptic  when  c < 0 , and 
(2.13)  with  initial  conditions  actually  is  ill-posed. 
However,  this  is  due  to  the  fact  that  we  are  considering 
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long  wave  perturbations.  A similar  result  occurs  in 
Whitham's  analysis  of  the  Stokes  wave  problem  [12]:  he 
finds  that  the  slow  modulation  equations  are  elliptic. 

In  the  one-dimensional  problem  this  instability  leads 
to  solitons.  Here  we  expect  to  find  lump-type  solutions. 

Conversely,  if  b > 0 , then  c > 0 , and  these  perturba- 
tions merely  propagate  without  change  along  the  crest  of  the 
soliton.  Thus,  the  stability  of  the  soliton  solution  of 
(1.8)  with  respect  to  transverse  perturbations  is  deter- 
mined by  exactly  the  same  criterion  that  determines  the 
(transverse)  stability  of  the  KdV  soliton  obeying  (1.7); 
i.e.,  whether  the  linearized  phase  speed  is  a maximum 
or  a minimum  at  k = 0 . 
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III.  DERIVATION  OF  THE  TWO-DIMENSIONAL  EQUATION 

In  this  section,  we  derive  (1.8)  in  the  context  of 
long  internal  waves.  For  the  sake  of  definiteness,  it 
is  convenient  (but  not  necessary)  to  think  in  terras  of 
a background  density  distribution  of  the  form 

p (z)  = p {1  + 6 exp (-z /h) } , z > 0 , (3.1) 

For  the  problem  of  waves  propagating  along  the  oceanic 
thermocline,  z = 0 denotes  a plane  of  symmetry  (where 
the  vertical  velocity  vanishes  for  the  modes  of  interest) 
at  the  center  of  thermocline,  h is  a measure  of  the 
thermocline  thickness,  the  ocean  surface  is  taken  at 
z = + °°  , and  typically  6 ~ 10" 2 . ’’Long"  disturbances 

have  horizontal  scales  much  larger  than  h , and  (1.8) 
results  if  they  are  also  much  smaller  than  the  distance 
to  the  surface.  Later,  we  will  indicate  in  context  what 
modifications  to  this  derivation  are  necessary  for 
density  distributions  that  differ  significantly  from  that 
in  (3.1). 

The  equations  of  motion  of  an  incompressible,  non- 
diffusive  fluid  are 

V • u = 0 , 

P Cl  = -Vp  - g . (3.2) 
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where  u = (u,v,w),  g = g(0,0,l)  and  D/Dt  = 9/9t  4-  u*V  . 
In  this  problem,  we  require 

w = 0 at  z = 0 , |u|  -*  0 as  z ■*  <*>  (3.3) 

In  the  absence  of  any  fluid  motion, 

u = 0 , 

P = PD(z)  , given  , 

P = P0(z  ) = -g  JT  Pcd;  . (3.  A) 

The  required  strategy  for  this  problem  can  be  seen 
by  looking  for  infinitesimal  perturbations  of  the  form 

<t>(z)eik(x~c(k)t) 


around  this  undisturbed  state.  The  result  is  an  equation 
for  the  vertical  velocity  of  the  form 

J?)  - {£  1ST  <z>  + kVz>}  * ’ 0 O-  5) 
For  realistic  density  profiles,  a natural  length  scale  is 


h = 


P0(z)  - Pc 


~S£ 


(z) 


min 


Hence,  the  waves  in  (3. A)  are  long  if 


| k | h <<  1 


(3.6) 


(3.7) 


In  what  follows,  we  shall  assume  that  |p^(z)|  is  mono- 
tonically  decreasing,  and  is  finite  as  z -*■  0 . Assuming 
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c2(k)  remains  finite  as  k -*■  0 , it  follows  that  there 
is  an  "inner  region"  where  (z/h)  is  the  appropriate 
variable  and  the  last  term  in  (3.5)  can  be  neglected, 
to  leading  order.  In  this  region,  the  solution  can  be 
expanded  as 

~ 4> o (z ) + |k|<t>^(z)  , 

c 2 - c 2 + I k I c£  (3.8) 

On  the  other  hand,  for  fixed,  small  k and  sufficiently 
large  z , there  is  an  "outer  region"  where  density  varia- 
tions can  be  neglected  and  (3.5)  reduces  to 

4>"  - k2 4>  ~ 0 . 


The  appropriate  variable  in  this  region  is  |k|z(=cz/h)  , 
and  the  solution  of  interest  is 

4>  - Ae"  I k I z . (3.9) 


The  region  of  overlap  is  found  by  solving  for  z from 
the  relation: 


1 dpo  ..  k2c2 (k) 
Pn(z)  ‘cTF  g 


(3.10) 


For  the  density  distribution  given  by  (3.1),  this  occurs 
as  z/h  -*•<*>  in  the  inner  region,  and  as  |k|z  -*■  0 in  the 
outer  region.  For  other  types  of  smooth  density  distribu- 
tions, the  matching  can  occur  in  a neighborhood  of  a 
finite  z/h  , but  always  as  |k|z  -*•  0 . Provided  dpQ/dz 
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is  continuous,  <f>  and  its  derivative  must  match  in  this 
region  of  overlap,  and  it  follows  from  the  structure  of 
(3.8)  and  (3.9)  that  the  required  conditions  are 

d6  d$, 

♦o  " A ■ nr  " 0 • - -A  • (3-n) 

These  matching  conditions  provide  the  additional  information 
required  to  specify  the  solutions  4>q  , <p  ^ completely. 

If  dpQ/dz  is  discontinuous,  there  may  be  no  region 
of  overlap,  and  then  the  matching  must  be  done  along  the 
isopycnic  surface  (surface  of  constant  density)  where  the  dis- 
continuity occurs.  The  required  matching  conditions  now 
are  that  the  normal  velocity  and  pressure  are  continuous 
across  this  surface.  In  the  linearized  problem,  these 
reduce  to  matching  <t>  and  d4>/dz  , but  not  in  the  non- 
linear problem.  We  will  discuss  this  point  in  more  detail 
after  first  deriving  (1.8)  for  a density  distribution  like 
that  in  (3.1),  where  dpQ/dz  is  continuous. 

Returning  now  to  (3.2)  and  (3.3),  we  apply  these 
ideas  to  the  nonlinear  problem.  Small  amplitude  waves 
which  are  long  in  the  sense  described  above,  traveling 
in  the  x-direction,  are  obtained  by  using  the  following 
scaling  in  the  lower  layer  ( z/h  - 0(1))  : 
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£ = e(x  - cQt) /h 
n = e3/2(y/h) 

4 = z/h 

t = e2cot/h 

u - ccou1(^,n,C,T)  + e 2 C0U2  . 
V ~ e 3 ^ 2 C0V1  (£ , h » C , x) 

w ~ e2cow1(5,n,C,T)  + £3cow2  • 
p ~ PD(C ) + ep00P1(C,n,C,T) 

P - P0(O  + (cpooc2)p1(^,n,i;(T) 


) (3.12) 


These  scales  were  used  by  Ono  [11],  except  for  y , which 
is  based  on  (1.9),  and  v , which  is  required  for  momentum 
balance.  To  leading  order,  after  using  (3.4),  (3.2)  yields 


9pl  dpo 

tt  + wi  ar  _ 0 

3u,  9w, 

TT  + TT  = ° 


9u 

(po(^>/pco)  H 


9px 

w~ 

(p«(0/p  J 


1 


9p. 

W 


= 0 

= 0 


9v  -I  3Pi 

"3 — = ® 
3^  9n 


(3.13) 


We  define 


«i  - - ||  (?,  n,T)4>(0 


(3.14) 


and  for  future  purposes,  we  require  that  9f/9£  -*■  0 as 
|£|  ■*  <*>  . Then  (3.13)  provides  the  equation  that  deter- 
mines the  vertical  structure  of  the  waves: 
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71  {po'c)  3s}  ■ f|  " 0 ■ <315> 

It  also  determines  the  other  unknowns  in  terms  of  f and  <t>  . 

As  *;-*■»,  P0  "*■  Poo  and  0 approaches  a linear 
function  of  c ; (i.e.,  <J>"  -*■  0)  . The  only  bounded  solution, 
fore,  satisfies  the  boundary  conditions: 


0(0)  = 0 , 


^ ■+■  0 as 


(3.16) 


in  accord  with  (3.11) 


an  eigenvalue  problem  for  (1/c^) 


Now  (3.15)  and  (3.16)  constitute 
2 ' There  are  infinitely 
many  solutions  to  this  problem,  but  the  lowest  eigenvalue, 
whose  eigenfunction  does  not  vanish  in  (0,®)  , is 
ordinarily  the  solution  of  physical  interest.  We  normalize 
this  eigenfunction  by  choosing  0(°°)  = 1 . 

The  relevant  equations  at  the  next  order  are: 


9Po  dp  if  a r 

-p»  ar  + w2  dr  = - It  ♦ p;  + f k u'^po  * 


'N 


Su2  3Vj_  9w2 

7T~  9n  8; 


3u, 

'po<5)  ITT 


+ P. 


= 0 


3p 


2 _ 3f  n 

3£  “ 3t  Po* 


3f 


+ f n + pD< 

^p? 

+ (8h/co)p2  = 0 ’ 


- p. 


'} 


J 


and  w2*0  at  5*0 

These  equations  reduce  to  a single  nonhomogenous 
equation  for  w2  , the  homogenous  portion  of  which  is 


there- 


(3.16) 
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identical  to  that  in  (3.15).  The  equation  has  a solution 
only  if  the  nonhomogenous  terms  satisfy  an  orthogonality 
condition : 


p&(w2*'  - 4>w2 ) 


oo 

/ 


00 

/ 


(pj'y*  ^ * f2  M + 1^ i 


9f 


•/ 


+ I (3  (P0<f  • 2 ) • - 2po<p'<p"  - 2(p0«*hJ)")  ' >4>  d;  • f || 
u (3.18) 

Again  we  note  that  if  (3.10)  should  hold  at  a finite  ; , 


say  hQ  , then  "®"  in  (3.18)  should  be  replaced  by 


In  the  upper  layer,  the  fluid  is  nearly  homogenous 
and  it  is  assumed  that  the  only  motion  is  that  excited 
by  waves  in  the  lower  layer.  Therefore,  the  equations 
of  motion  must  reduce  to  LaPlace's  equation  to  leading 
order.  Moreover,  the  horizontal  scales  should  be  chosen 
to  be  consistent  with  those  in  the  lower  layer,  but  the 
vertical  scale  changes.  Hence,  in  this  region,  we  define 


Z=ez/h’-eC  , 
w = e 2cqW  (C ,n ,Z  , t , e) 


(3.19) 


and  the  scaling  of  the  other  independent  variables  is 
taken  from  (3.12).  The  equation  for  the  vertical  velocity 
Z > 0)  is 


9 2W  9 2W 

jr7  + bzt 


0(c) 


(3.20) 
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subject  to  the  conditions  that 

W -*■  0 as  Z -♦  00  , 

W -*■  0 as  | £ | -*■  00  , 

W=-|i-(C,n,T,e)  on  Z = 0 

This  last  condition  is  derived  from  the  requirement  that 
the  two  representations  of  the  vertical  velocity  must 
match  (as  functions)  in  the  region  of  overlap,  which  in 
this  case  is 

ep  z/h  = 0(1)  , 0 < p < 1 . (3.22) 


By  construction , the  two  representations  of  the 
vertical  velocity  match.  Matching  the  vertical  derivatives 
of  these  functions  as  well  insures  the  unique  analytic 
continuation  of  W , as  required.  From  (3.16)  and  (3.24), 
this  implies 
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£ -+co 


(3.25) 


The  last  step  is  to  substitute  (3.25)  into  (3.18)  and 
integrate  by  parts  appropriately,  using  the  fact  that 

4>"  0 as  c “ (3.26) 


The  result  has  the  form  of  (1.8),  where  u is  identified 
with  f , 


b = 


P04>’2d; 


(3.27) 


and  cQ  = 1 because  of  our  scaling  in  (3.12).  In  this 
application  b > 0 , and  no  disturbance  can  arrive  at 
£ * before  the  solution  of  (1.8).  Hence,  it  is  ap- 
propriate to  take  xq  = +°°  in  (1.8),  so  that  the  integral 
term  becomes  / f d£ ' . It  follows  from  (3.12),  (3.13) 

on 

and  (3.14)  that  f , the  solution  of  (1.8),  represents  both 
the  horizontal  fluid  velocity  and  the  vertical  displacement 
of  an  isopycnic  surface.  Moreover  the  soliton,  (1.6),  is  not 
unstable  with  respect  to  long  transverse  perturbations,  as 
shown  in  Section  II. 

The  integrals  in  (3.23)  and  (3.24)  coverge  only  if 
-*■0  as  U|  - o°  . With  xq  = -H>°  , there  is  no  apparent 
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difficulty  in  requiring  that  the  solution  of  (1.8)  should 
vanish,  along  with  its  derivatives,  as  £-*•+“.  However, 
even  if  f (i.e.,  u)  and  all  of  its  derivatives  vanish 
initially  as  £ -«  , it  is  evident  from  (1.8)  that  f 

will  not  remain  zero  there  unless 


00 


9 2 f 
9n2 


d£ 


= 0 


(3.28) 


This  additional  constraint  on  the  solution  can  be  inter- 
preted in  terms  of  the  vertical  component  of  vorticity, 


U) 


9u  9v 
9y  9x  ‘ 


(3.29) 


Using  (3.12),  (3.13)  and  (3.14),  one  can  show  that  both  terms 

in  (3.29)  are  of  the  same  order  of  magnitude  (in  e ) , and  that 

9f 

both  are  proportional  to  . Thus,  (3.28)  is  satisfied 

if  the  total  (i.e.  , integrated  in  £ ) <d- component  of  vorticity 
is  constant  in  n . 

We  close  this  section  by  outlining  the  analysis  required 
for  discontinuous  density  profiles.  The  main  point  is  that 
when  either  pQ  or  p^(z)  is  discontinuous,  there  is  no 
matching  region,  but  rather  a sharp  interface.  In  these 
cases,  we  must  match  normal  velocity  and  pressure  along  the 
interface  which,  in  turn,  has  to  be  found.  The  essential 
steps  are  as  follows.  We  define  the  interface  z « £(x,y) 
by 


Dz 

Bt 


u- V; 


W 


(3.30) 
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If  we  call 

A 

F = z - c = z - 1 - d1  - e2C2-  • • , 

then  using  the  unit  normal,  n = VF/|VF|  , the  problem 
becomes  specified  by  the  matching  conditions 


u • ft 


lower 


= u*n 

upper 


(3.31) 


Pressure 


lower 


Pressure 


upper 


(3.32) 


We  then  successively  (in  powers  of  c ) satisfy  in  order: 
(3.2  1),  (3.30),  (3.32),  from  which  we  obtain  unique  func- 
tions for  the  velocity,  interface,  and  pressure.  In  all 
cases,  we  find  (1.8)  with  a,b  given  by  (3.27)  (and  °° 
replaced  by  h = 1 ) . In  the  case  where  we  have  two 
homogenous  fluids  of  differing  densities,  a = 3/2 

b = p /2p  where  p , p are  the  respective  densities 

00  o 00  o 

of  the  upper  and  lower  fluids. 

When  pQ(z)  is  continuous,  but  Pq(z)  is  discontin- 
uous, then  the  effect  of  the  interface  is  weaker.  To 
obtain  (1.8),  we  need  only  use  leading  order  results 
from  (3.31);  i.e.,  continuity  of  vertical  velocity  (to 
higher  order,  however,  (3.3l)  itself  must  be  employed). 

If  Py(z)  is  also  continuous  then  a matching  region  exists, 
and  the  analysis  presented  at  the  beginning  of  this  sec- 
tion applies. 
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Finally,  it  should  be  noted  that  the  wave  modes  of 
interest  here,  which  leave  the  center  of  the  thermocline 
undisturbed,  travel  slower  than  either  the  surface  wave  or 
the  fastest  internal  wave  with  the  same  horizontal  wave- 
number.  The  analysis  presented  here  assumes  that  if  there 
is  any  energy  in  these  other  wave  modes,  the  relevant  group 
velocities  greatly  exceed  cq  , so  that  there  is  no  coupling 
between  the  long  internal  waves  considered  here  and  packets 
of  shorter  waves  in  these  other  modes.  If  there  is  significant 
coupling,  the  equations  are  more  complicated  than  (1.4)  and 
(1.8). 
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ABSTRACT 


We  consider  the  evolution  of  packets  of  water  waves 
that  travel  predominantly  in  one  direction,  but  in  which 
the  wave  amplitudes  are  modulated  slowly  in  both  horizontal 
directions.  Two  separate  models  are  discussed,  depending 
on  whether  or  not  the  waves  are  long  in  comparison  with 
the  fluid  depth.  These  models  are  two-dimensional  generaliza 
tions  of  the  Korteweg-deVries  equation  (for  long  waves) 
and  the  cubic  nonlinear  Schrodinger  equation  (for  short 
waves).  In  either  case,  we  find  that  the  two-dimensional 
evolution  of  the  wave  packets  depends  fundamentally  on  the 
dimensionless  surface  tension  and  fluid  depth.  In  particular 
for  the  long  waves,  one-dimensional  (KdV)  solitons  become 
unstable  with  respect  to  even  longer  transverse  perturba- 
tions when  the  surface  tension  parameter  becomes  large 
enough,  i.e.,  in  very  thin  sheets  of  water.  Two-dimensional 
long  waves  ("lumps")  that  decay  algebraically  in  all 
horizontal  directions  and  interact  like  solitons  exist 
only  when  the  one-dimensional  solitons  are  found  to  be 
unstable . 

The  most  dramatic  consequence  of  surface  tension 
and  depth,  however,  occurs  for  capillary-type  waves  in  suf- 
ficiently deep  water.  Here  a packet  of  waves  that  are 
everywhere  small  (but  not  infinitesimal)  and  modulated  in 
both  horizontal  dimensions  can  "focus"  in  a finite  time, 
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producing  a region  in  which  the  wave  amplitudes  are  finite. 
This  nonlinear  instability  should  be  stronger  and  more 
apparent  than  the  linear  instabilities  examined  to  date; 
it  should  be  readily  observable. 

Another  feature  of  the  evolution  of  short  wave  nackets 
in  two -dimens ions  is  that  all  one-dimensional  solitons  are 
unstable  with  respect  to  long  transverse  perturbations. 
Finally,  we  identify  some  exact  similarity  solutions  to  the 
evolution  equations. 


ill 


1. 


INTRODUCTION 
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Our  understanding  of  the  evolution  of  surface  water 
waves  of  moderate  amplitude  has  increased  significantly 
within  the  last  decade  or  so.  The  evolution  in  one  spatial 
dimension  of  a packet  of  inviscid  waves  of  sufficiently  small 
amplitude  is  governed  by  linear  equations  on  a short  time 
scale,  and  by  either  the  Korteweg-deVries  (KdV)  equation 


u.  + uu  + 

t X 


u 


XXX 


0 


or  the  cubic  nonlinear  Schrddinger  equation 


iAt  + Axx  + o | A | 2 A = 0 


(1.1) 


(1.2) 


on  longer  time  scales,  depending  on  whether  or  not  the 
typical  wavelengths  are  large  in  comparison  with  the  fluid 
depth.  In  (1.2)  and  throughout  this  paper,  o = ±1  , and 
represents  an  irreducible  choice  of  signs.  Both  of  these 
equations  can  be  solved  exactly  as  initial  value  problems, 
using  inverse  scattering  transforms  (1ST;  an  account  of  1ST 
can  be  found  in  [1]).  In  situations  in  which  viscous  ef- 
fects are  felt  on  an  even  longer  time  scale,  these  theories 
(or  viscously-corrected  versions  of  them)  predict  with  very 
reasonable  accuracy  the  evolution  of  waves  over  quite  long 
distances  in  wave  tanks  (Hammack  and  Segur , [2],  [3];  Yuen 
and  Lake , [ A]) . 
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Outside  of  specially  designed  tanks,  surface  waves 
ordinarily  evolve  in  two  spatial  dimensions  and  here  the 
theory  is  much  less  complete.  A two-dimensional  generaliza- 
tion of  (1.1)  for  nearly  one-dimensional  long  waves  was 
given  by  Kadomtsev  and  Petviashvili  [5]  in  the  form: 


u_  + uu  + ou 
t x xxx 


oo 

/ 


X 


u dx  = 0 

yy 


(1.3) 


Results  by  several  authors  indicate  that  (1.3)  is  of  IST- 
type,  but  a complete  method  of  inverse  scattering,  analogous 
to  that  in  one  spatial  dimension,  has  not  yet  been  developed. 

Two-dimensional  generalizations  of  (1.2)  were  derived 
by  Zakharov  [6],  Benney  and  Roskes  [7],  Davey  and  Stewartson 
[8],  and  Djordjevic  and  Redekopp  [9].  All  of  these  analyses 
followed  approximately  the  same  lines.  The  problem  was 
also  studied  by  Hayes  [10],  using  somewhat  different  methods. 
The  most  general  analysis  was  by  Djordjevic  and  Redekopp, 
who  included  the  effects  of  gravity,  surface  tension  and 
arbitrary  depth  to  get  a system  that  can  be  reduced  to 


1At  + °lAxx  + Ayy  ■ °2IAI 2A  + V ■ 


a$  +<j>  = -b ( I A I 2 ) 

xx  yy  1 1 x 


(1.4) 


where  (a.b.o^.c^)  depend  on  the  (dimensionless)  fluid  depth 
and  surface  tension.  In  the  long  wave  limit,  (1.4)  reduces 
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to  one  of  the  problems  that  Ablowitz  and  Haberman  [11]  had  shown 
were  of  IST-type.  As  with  (1.3),  beyond  identifying  the  ap- 

I 

propriate  linear  scattering  problem  and  obtaining  special  solu- 
tions, no  general  inverse  scattering  theory  has  yet  been  developed. 

In  these  two  cases,  (1.3)  and  the  long  wave  limit  of  (1.4), 
one  can  reasonably  anticipate  that  the  necessary  inverse  scat- 
tering theory  eventually  will  be  developed,  and  that  the  general 
solutions  of  (1.3)  and  (1.4),  as  initial  value  problems,  will 
become  available.  In  these  cases,  the  two-dimensional  problem 
should  eventually  be  solved  to  the  extent  that  the  one-dimensional 
problem  is  now.  However,  as  discussed  in  §5,  we  conjecture  that 
(1.4)  cannot  be  solved  by  inverse  scattering  transforms  over 
the  entire  range  of  parameters  and  that  the  general  two-dimen- 
sional problem  cannot  be  solved  in  a manner  analogous  to  that 
in  one  dimension. 

The  purpose  of  this  paper  is  to  identify  some  important  re- 
sults regarding  (1.3)  and  (1.4),  and  to  suggest  the  role  that 
they  play  in  the  solution  of  initial  value  problems.  A major 
result  of  this  study  is  the  dramatic  effect  that  surface  tension 
can  have  upon  the  dynamics  of  the  wave  motion.  A summary  of 
these  results,  and  an  outline  of  the  paper  is  as  follows. 

(§2)  The  derivations  of  (1.3)  and  (1.4)  from  the  physical 
problem  of  water  waves  are  discussed.  These  equations  are  well 
established  in  the  literature,  but  the  question  of  what  boundary 
conditions  and  other  constraints  are  required  to  make  the  prob- 
lems well-posed  is  still  open.  We  show  that  the  original  prob- 
lem selects  certain  side  conditions  as  "natural." 

^ Which  conditions  are  appropriate  depends  on  the  dimensionless 

I 
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surface  tension  and  depth.  In  this  section  we  also  consider 
the  physical  interpretation  of  an  infinite  set  of  conserva- 
tion laws . 

(§3)  The  role  that  one-dimensional  soliton  solutions 
can  play  in  the  two-dimensional  problems  is  examined  (i.e., 
stability  of  solitons) . KdV  solitons  are  unstable  in  (1.3) 
when  a = -1  , which  occurs  in  sufficiently  thin  sheets  of  • 
water  (i.e.,  large  enough  surface  surface  tension  coefficient). 
For  zero  surface  tension  a = +1  , and  the  argument  does 
not  yield  instability.  When  solitons  are  unstable,  they 
cannot  be  viewed  as  the  asymptotic  (t  -*■  °°)  states  towards 
which  the  solution  evolves,  as  they  are  in  the  one-dimensional 
problem.  In  this  case,  "lump"  solutions  exist  and  may  play 
an  asymptotic  role  analogous  to  that  of  one-dimensional 
solitons . 

Zakharov  and  Rubenchik  [12]  showed  that  for  the  one- 
dimensional cubic  nonlinear  Schrodinger  equation  all  one- 
dimensional solitons  are  unstable.  These  results  apply  to 
the  deep  water  limit  of  (1.4).  We  extend  their  analysis  to 
demonstrate  the  equivalent  results  in  the  case  of  finite 
depth. 

(§4)  The  most  dramatic  effect  of  strong  surface  ten- 
sion is  focusing . A wave  that  is  large  enough  (in  a certain 
integral  sense)  focuses  at  a particular  point  in  space 
after  a finite  time.  Here  there  is  no  asymptotic  (t  -*■  «°) 
state,  because  the  solution  of  (1.4)  develops  a singularity 
in  a finite  time.  Focusing  provides  a mechanism  by  which 
a field  of  relatively  small  amplitude  waves  produces  a local 
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region  in  which  the  amplitudes  are  large.  Focusing  is  a 
potentially  important  mechanism  in  the  redistribution  of 
energy  within  the  spectrum;  it  should  be  readily  measurable. 

(§5)  We  consider  the  question  of  the  complete  integra- 
bility  of  (1.4).  Moreover,  we  exhibit  some  special  solutions 
that  are  not  one-dimensional,  and  are  candidates  for  asymp- 
totic states  in  the  two-dimensional  problem. 


2 . RELEVANT  EVOLUTION  EQUATIONS 

The  classical  problem  of  water  waves  is  to  find  the 
the  irrotational  motion  of  an  inviscid,  incompressible, 
homogeneous  fluid,  subject  to  the  forces  of  gravity  and 
surface  tension.  The  fluid  rests  on  a horizontal  and  im- 
permeable bed  of  infinite  extent  at  z = -h  (h  may  be 
finite  or  infinite),  and  has  a free  surface  at  z = c(x,y,t)  . 
The  fluid  has  a velocity  potential,  <f>  , which  satisfies 

V24>  — 0 -h  < z < £(x,y,t)  (2.1) 

It  is  subject  to  boundary  conditions  on  the  bottom,  z = -h  : 

<PZ  = 0 ; (2.2) 

and  along  the  free  surface,  z = C : 


DC  _ 
Dt 


ktt  = C-  + = 


x x 


y y 


i , ^v(1+0  + cvv(i+0  - 2^xv^v 

gc  + (fv  + ^-1  V4>  1 2 = T 2 TL VLJLJL 


(2.3) 


d + ^ + 9 


Here  g is  the  gravitational  acceleration , and  T is  the 
ratio  of  surface  tension  coefficient  to  fluid  density.  We 
note  that  the  linearized  dispersion  relation  for  this  sys- 
tem is 


oj2  = (g<  + K3T)tanh  xh 


(2.4) 


In  two  dimensions,  one  should  interpret  k = A 2 + £2  in 
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2 . 1 KdV  Limit 

The  solution  of  (1.3)  provides  an  approximate  solution 
to  these  equations  that  is  valid  when  the  initial  disturbnace 
consists  primarily  of  nearly  one-dimensional  long  waves  of 
small  amplitude.  To  be  precise,  let  ic  = (k,£.)  be  the 
horizontal  wavenumber  characteristic  of  the  disturbance. 
Orient  the  horizontal  coordinate  system  such  that  the  x- 
direction  is  the  principal  direction  of  wave  propagation. 

Let  a denote  the  characteristic  amplitude  of  the  distur- 
bance. Then  we  need: 

(i)  small  amplitudes, 


e = a/h  <<  1 ; 


(2.5a) 


(ii)  long  waves, 


(<h)2  <<  1 ; 


(2.5b) 


(iii)  nearly  one-dimemensional  waves, 


(£/«) 2 <<  1 . 


C2.5c) 


The  KdV  equation  (1.1)  results  when  the  first  two  effects 
balance  in  truly  one-dimensional  problems,  and  (1.3)  results 
when  all  three  effects  balance: 


(<h)2  = 0(e)  , (2 . 5d) 

(*/k)2  = 0(e)  . (2 . 5e) 


I 

i 
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Under  the  assumptions  of  (2.5)  a first  approximation 
of  (2.1)  - (2.3)  reduces  to 

- gh  0 - 0(e)  (2.6) 

Thus,  to  lowest  order,  the  solution  of  (2.1)  - (2.3)  may 
be  approximated  by 


C ~ eh[f1(x  - /gfit ; y)  + t'2(x  + /ght;y)]  , (2.7) 


where  f^,f2  are  known  in  terms  of  the  initial  dada. 
Throughout  this  paper,  we  are  interested  in  problems  where 
the  initial  disturbances  are  localized,  and  it  is  then  con- 
venient to  assume  a fortiori  that  the  physical  quantities 
have  compact  support  initially.  In  this  case,  it  is  easy 
to  show  the  f^  and  f2  in  (2.7)  have  compact  support  as 
well . 

To  go  to  higher  order,  we  define  scaled,  dimension- 
less variables: 


r - /c(x-/ght)/h 
n - ey/h 


s = /F(x+/gHt)/h 
t = e/ght/h 


u - fj_  , v = f2 

T = T/gh?  . 


1 


(2.8) 


J 


Now  we  look  for  solutions  of  the  form  r,  - eh[u(r,T,n)  + v(s,i,n)] 
i.e.,  we  use  the  method  of  multiple  scales.  To  eliminate 
secular  terms  at  the  next  order,  we  find 
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(2u,  + 3uu  + (*•  - f)u  ) + u = 0 , 

i r 3 rrr  rm 

r 

(2v  - 3w  - (i  - T)v  ) - v = 0 

r s 3 ' sss7  n n 

s 

The  equation  given  by  Kadomtsev  and  Petviashvili  [5]  is 
in  this  form. 

For  most  circumstances  of  interest  in  water  waves, 


(2.9) 


J - T > 0 


(2.10) 


and  it  follows  from  (2. A)  that  the  linearized  phase  speed 
is  a (local)  maximum  at  k = 0 . Thus,  the  waves  governed 
by  (2.9)  travel  faster  than  their  neighbors  (in  tc-space) 
and  there  should  be  no  disturbance  as  r -*•  , or  s -*•  -«• 

Consequently,  (2.9a)  may  be  integrated  to 

oo 

2u  + 3uu  + (i  - T)u  - I u dz  *=  0 , (2.11) 

t r 3 ' rrr  I nn  * 

r 

with  a similar  equation  for  (2.9b).  This  is  now  in  the 
form  of  an  evolution  equation  for  u , as  is  (1.3).  For 

A 

very  thin  sheets  of  water  (i.e.,  T large  enough) 

(2.10)  is  false,  the  long  waves  travel  slower  than  their 
neighbors,  and  the  integral  in  (2.11)  should  be  over  (-°°,r)  . 

Given  (2.10),  there  is  no  apparent  difficulty  in 
requiring  that  u should  vanish,  along  with  its  derivatives, 
as  r -*•  +00  . However,  even  if  u and  all  of  its  derivatives 
vanish  initially  as  r -*•-<»,  it  is  evident  from  (2.11) 
that  u will  not  remain  zero  there  unless 
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/u  dr 

nn 


0 


— oo 


(2.12) 


Since  u is  the  derivative  of  a velocity  potential,  (2.12) 
is  automatically  satisfied  at  the  initial  instant.  Indeed, 
for  the  linearized  form  of  (2.11),  (2.12)  is  a constant 
of  the  motion,  and  it  is  sufficient  to  require  it  initial- 
ly. 

The  constraint  in  (2.12)  has  a simple  physical  inter- 
pretation. One  can  identify  j u(r,n,T>dr  as  the  total 
mass  of  the  wave  in  a thin  strip  at  n . Then  (2.12)  as- 
sures that  the  transverse  derivative  of  mass  is  constant, 
and  this  prevents  a net  flow  of  mass  to  (or  from)  any 
particular  strip. 

There  are  several  indications  that  (2.11),  or  (1.3),  is 
of  IST-type.  Dryuma  [13]  has  identified  an  appropriate 
linear  scattering  problem  for  (1.3);  Zakharov  and  Shabat 
[14]  have  related  special  solutions  to  a linear  integral 
equation;  Chen  [15]  found  a Backlund  transformation; 

Satsuma  [16]  has  obtained  "N  soliton,”  but  nonlocalized , 
solutions  by  direct  methods.  In  Section  III  we  discuss 
localized  lump  solutions.  However,  as  mentioned  earlier, 
no  complete  1ST  method  has  been  developed  for  (1.3)  to 
date. 
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2 . 2 The  Nonlinear  Schrodlnger  Limit 

Let  us  now  consider  the  derivation  of  (1.4)  from 
(2.1)  - (2.3).  Here  we  are  following  a packet  of  nearly 
one-dimensional  waves,  traveling  in  the  x-direction, 
with  an  identifiable  (mean)  wavenumber,  k = (k,£)  . 

We  denote  the  maximum  variation  in  k within  the  packet 
by  6k  . To  derive  (1.4)  we  need: 

(i)  small  amplitudes, 


e = Ka  1 

(ii)  slowly- varying  modulations, 

6k/ k <<  1 


(2.15a) 


(2.15b) 


(iii)  nearly  one-dimensional  waves, 

|£|/k  <<  1 (2.15c) 

(iv)  balance  of  all  three  effects, 

6k/<  = 0(c)  (2 . 15d) 

\1\/k  - 0(c)  (2 . 15e) 


The  dimensionless  depth,  kh  , can  be  finite  or  infinite, 
but  to  avoid  the  shallow  water  limit  (and  KdV) , we  need 


(kh) 2 >>  e . 


(2.16) 
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In  this  limit,  the  solution  of  the  lowest  order  (linear) 
problem  is 


(< 


E - 


V 


cosh  k (z+h) 
cosh  kh 


[A  exp(i0)  + (*) ] + const. 


(2.17a) 


where  (*)  denotes  complex  conjugate, 

0 = kx  - u>(k)t  , (2.17b) 

and  u>  (k)  is  given  by  (2.4).  To  go  to  higher  order,  we 
introduce  slow  (dimensional)  variables  (again,  using  the 
method  of  multiple  scales) , 

xj  = cx  , yi  = ey  , t,  = ct  , t2  = e2t  , (2.18) 


and  expand  4)  and  c,  : 


4>  ~ e 


]*(xi,yi,  ti.tz)  + C°coS{iThh)'  [A(xi  ,yi  ,tj  ,t2)exp(ie) 


+ (*)]  + 0 (e 2 ) 


(2.19) 


4 = c < Cnexp(iO)+(*) 


+ 0 (e 2 ) , Ci  i - 


10) 


In  order  to  derive  (1.4),  these  expansions  must  be  carried 
out  to  0(e3)  . The  variations  allowed  in  A reflect 
the  fact  that  this  is  a wave  packet,  rather  than  a uniform 
wavetrain,  and  Q provides  the  mean  drift  current  generated 
by  the  packet.  In  what  follows  we  shall  only  discuss  the 
secular  effects  that  the  higher  order  terms  have  on  <t>  , 
and  A ; details  can  be  found  in  [7]- [9]. 
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At  the  next  order  of  approximation,  a secular  condi- 
tion requires  that  the  wave  packet  travel  with  its  linear 
group  velocity, 


+ C (k)  - 0 

otj  gv  tod 


(2.20) 


where  C * dw/d^  . On  this  same  time  scale,  i>  satisfies 
g 

a forced  wave  equation, 


)H  ,1V*  , 
ft!  " gh(jxl  9yT J 


k.B,  |A|! 


(2.21) 


where 


IVV 

6i  = — sech2kh  + 2/(l+T)  , 


T = k2T/g  = (kh) 2T 


The  solution  of  (2.21)  changes  dramatically,  depending  on 


whether  or  not 


gh  > . c: 

If  the  ratio  C / /gTi  is  interpreted  as  the  "Mach  number" 
g 

of  the  wave  packet,  then  (2.22)  is  the  condition  for  "sub- 
sonic" flow.  In  this  case,  if  A has  compact  support, 
then  4>  has  a forced  component  that  travels  with  speed 

C (i.e.,  it  satisfies  (2.20)),  and  a free  component 
g 

that  radiates  outward  with  speed  /gK  , and  is  O(tl^)  as 
tj  -*•  o°  . Hence  with  (2.22),  as  ti  -*•  00  , wp  find  that  * 


(2.22) 


satisfies  both  (2.20)  and 


3 2 4>  , 3 2 * 

a 3x1  + 9yT 


gH  6l  3x7  I*'’ 


(2.23) 
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where 

a = (gh  - Cp/gh 


along  with  the  boundary  condition  that  4-  vanishes  as 
(xi  + y2)  -*■  00  . These  are  the  boundary  conditions  pre- 

scribed by  Davey  and  Stewartson  [8],  and  they  are  correct 
without  surface  tension. 

If  the  effects  of  surface  tension  are  strong  enough, 
(2.22)  fails  and  the  flow  is  "supersonic."  Now  even  if 
A has  compact  support,  4 and  its  derivatives  are  non- 
zero along  "Mach  lines"  that  emanate  from  the  support  of  A . 
In  the  limit  ti  00  , 4 satisfies  both  (2.20)  and  (2.23) 
as  before.  However,  the  appropriate  boundary  conditions 
for  (2.23)  now  are  that  4 and  its  derivatives  vanish 
ahead  of  the  support  of  A (e.g.,  as  xi  -*•  °°)  , and  no  con- 
ditions as  Xj  . Hence,  in  general,  we  can  not  expect 

that  global  integrals  involving  4 will  converge. 

The  limit  ti  ■+■  00  is  of  interest  because  (1.4)  ap- 
pears when  one  eliminates  secular  terms  on  the  next  time 
scale,  t = 0(e“2)  . Carrying  this  out,  and  putting  the 
result  in  dimensionless  form,  we  define 


i = ek(x  - Cgt) 
t = e2 (gk)^t 
A = k2 (gk)"^A 


n = cky 


4 = k2(gk)"^4  1 


(2.24) 


and  find  that  A and  4 satisfy 


In  the  above  formulae,  all  functions  are  evaluated 
at  £*0  , since  we  are  considering  our  underlying  wavetrain 
to  be  propagating  purely  in  the  x-direction.  It  should  be 
noted  that  (2.25)  can  be  easily  scaled  to  (1.4)  where 
Oi  - sign  \ , o 2 m sign  x . a - a y/A2  and  b - 3PX 1 / x 2 I X I 


in  (1.4). 
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(2.25)  - (2.26)  are  equivalent  to  those  of  Djordjevic  and 
Redekopp  ([9],  their  equations  (2.12)  - (2.13))  except  for 
the  correction  of  a misprint.  If  the  initial  wave  packet 
is  local,  it  is  appropriate  to  require  that  A vanishes 
as  £2  + n2  -*•  00  . As  discussed  above,  the  appropriate 
boundary  conditions  for  <J>  depend  on  the  sign  of  a 
In  the  deep  water  limit,  (2.25)  reduces  to 

iA  + X Arr  + y A = x | A | 2 A , (2.27) 

T OO  ££  OO  n Aoo  I 1 ’ v ' 


where 


X 


OO 


% ( 1-6T-3T2  \ 
" "Boj  V~ T+T ) 


Moo  = 7^  <1  + 3^ 

8+T+2T2 

Xoo  = 7^  (l-2T)(l+f) 


The  appropriate  boundary  conditions  for  localized  initial 
data  are  that  A vanish  as  + n2  -*•  00 

The  character  of  the  solution  of  (2.25)  depends  funda- 
mentally on  the  signs  of  the  coefficients  in  the  equations. 
Figure  1 is  a map  of  parameter  space,  showing  where  these 
signs  change.  The  figure  is  that  of  Djordjevic  and  Redekopp 
[9]  , who  used  it  to  explain  the  various  regions  of 
stability/ instability  of  the  Stokes  wave.  Each  boundary 
line  corresponds  to  a simple  zero  of  a coefficient,  as 
shown,  except  for  the  two  curves  bounding  region  F . These 
two  curves  denote  singularities  of  v . In  a neighborhood 
of  each  of  these  two  curves,  phenomena  occur  on  a shorter  time 
scale  than  the  0(c"2)  scale  required  elsewhere;  cf.  [9], 
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If  we  take  the  long  wave  limit,  kh  -*•  0 , of  (2.25)  but 
keeping  c <<  (kh)2  , we  find  equations  which  are  of  IST- 
type.  We  discuss  this  further  in  Section  5.  Alternatively, 
the  long  wave  limit  in  which  c « 0((kh)2)  , where  (2.11) 
applies,  corresponds  to  the  lower  left-hand  corner  of  this 
figure  (kh  -*■  0 , T -*  0 , f fixed).  The  only  parameter  that 
changes  sign  in  this  limit  is  (1/3  - T)  , which  is  positive 
in  Region  A,  and  negative  in  Region  B.  The  uniformity  of  the 
limits  kh  -*■  0 , e -*■  0 has  been  discussed  in  [17,18,19]. 


2 . 3 Conservation  Laws 

Our  final  objective  in  this  section  is  to  give  a 
simple  physical  interpretation  for  an  infinite  set  of  con- 
servation laws.  It  is  well-known  that  the  equations  of 
water  waves  conserve  mass,  horizontal  momentum  and  energy. 

If  we  interpret  "mass"  as  the  mass  associated  with  the 
wave,  etc.,  then  these  conserved  quantities  may  be  repre- 
sented as  integrals.  In  one  dimension  (which  is  sufficient 
for  the  purpose  of  this  discussion)  we  have: 

Hass 

M - p j tdx  ; (2.28) 


Moment um 


(2.29) 
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Energy 

K.E.  = jj  j | V<J<  | 2 dz  dx 
-h 

P.E.  = % g J C*dx  (2.30) 

E = K.E.  + P.E. 

On  the  other  hand,  problems  that  have  been  solved 
exactly  by  1ST  possess  an  infinite  set  of  conservation  laws. 

For  example,  the  first  few  quantities  conserved  by  (1.2) 
are 

zl  m f lAl 2dx 

t2  = f (A*Ax  - A*A) dx  (2.31) 

Z3  = / ( I Ax I 2 ‘ Z lAl">dx 

There  has  been  some  speculation  about  the  proper  physical 
interpretation  of  this  infinite  set  of  conserved  quantities. 

We  offer  here  a very  simple  explanation.  We  have  seen  that 
(1.1)  - (1.4)  all  are  obtained  via  expansions  in  wave  am- 
plitude, c . From  this  viewpoint,  one  might  also  expand 
(for  example)  the  expression  for  the  mass  of  the  wave  in 
powers  of  e , to  obtain  a series  of  the  form 

oo 

M = p £ enCn  (2.32) 

Because  M is  constant  in  time,  it  follows  that  each 

coefficient,  C , is  also  constant, 
n 
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Because  one  generates  the  complete  series  for  4>  , C 
through  0(e3)  in  deriving  (1.4),  it  is  then  straight- 
forward to  compute  the  series  in  (2.32)  to  this  order.  In 
(2.19),  any  terms  involving  exp(ie)  can  be  shown  to  con- 
tribute only  at  higher  order,  using  integration  by 

parts  : 


ti 


. \ ie  , 

t2)e  dx 


dx 


dx 


This  process  can  be  repeated  as  many  times  as  Cm  can 
be  differentiated.  The  result  of  explicit  computation 
is 


M = eat Ii  + e2a2I2  + 0(e3)  ; 

m - e2b2I2  + 0(e3)  ; 

x (2.33) 

K.E.  = ec 1 1 1 + e 2c2 12  + 0(e3)  ; 

P.E.  - ec 1 1 1 + e 3c2I2  + 0(e3)  . 

I3  enters  at  0(e3)  . The  coefficients  (a^.b^.c^)  are 
unimportant  for  our  purpose.  The  momentum  starts  at  higher 
order  because  it  is  referred  to  a coordinate  system  travel- 
ing with  the  group  velocity  of  the  wave.  The  identity  of  the 
last  two  series  is  a statement  of  the  equipartition  of  the 
averaged  energy,  to  this  order.  It  is  not  true  that  I3  , I2 
and  I i represent  respectively  the  mass,  momentum  and  energy 
of  the  water  waves.  (Similarly,  the  first  three  conserved 
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quantities  for  KdV  are  not  respectively  the  leading  terms 
of  the  expansions  of  the  mass,  momentum  and  energy  of  the 
water  waves . ) 


I 
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3.  STABILITY  OF  SOLITONS 

The  primary  purpose  of  this  section  is  to  discuss  the 
stability  of  solitons  with  respect  to  transverse  perturba- 
tions . 

3.1  The  KdV  Limit 

Let  us  first  consider  the  long  wave  problem,  and  (2.11). 
The  one-dimensional  limit,  3/3 n = 0 , yields  KdV.  Here 
initial  data  on  compact  support  evolve  into  a finite  num- 
ber of  solitons,  ordered  by  amplitude,  followed  by  decaying 
oscillations  that  can  be  described  in  terms  of  a modulated 
similarity  solution.  The  decay  rate  of  the  oscillations 
is  not  uniform  in  space,  but  it  is  of  algebraic  order 
[20].  The  solitons  are  (theoretical) 
waves  of  permanent  form  when  separated  spatially  from 
other  waves.  They  represent  water  waves  that  decay  only 
due  to  viscous  effects.  A KdV  soliton  is  shown  in  Figure  2a. 
Both  the  solitons  and  the  decaying  oscillations  have  been 
observed  experimentally  [2,3]. 

Kadomtsev  and  Petviashvili  [5]  analyzed  the  stability 
of  a KdV  soliton  with  respect  to  long  transverse  perturba- 
tions in  (2.11).  They  found  that  the  soliton  is  unstable 
with  respect  to  such  perturbations  when  (2.10)  fails  (i.e., 
in  the  lower  left  corner  of  Region  B in  Figure  1) . The 
usual  situation  is  Region  A,  where  (2.10)  applies.  Here 
they  did  not  find  that  the  soliton  is  unstable. 


* 
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In  Region  B,  where  the  solitons  are  unstable,  the 
KdV  theory  is  of  limited  value.  Here  the  solitons  cannot 
represent  asymptotic  states,  as  they  do  in  the  one-dimen- 
sional problem.  Thus,  the  question  arises  as  to  whether 
(2.11)  has  any  other  special  solutions  that  might  act  as 
asymptotic  states  when  the  solitons  are  unstable.  The 
answer  is  not  known  definitively  at  this  time,  but  the  work 
it.  [21,22]  is  suggestive.  In  Region  B,  but  not  in 
Region  A,  (2.11)  possesses  "lump"  solutions.  Lumps  share 
many  of  the  important  properties  of  solitons: 

(i)  Each  is  a permanent  wave  whose  speed,  relative  to 
the  linearized  speed,  /gK  , can  be  made  proportional  to  its 
amplitude . 

(ii)  Solitons  are  localized  waves,  with  exponential  tails 
in  one  dimension;  lumps  are  localized  waves,  with  algebraic 
tails,  in  two  dimensions. 

(r'.ii)  Two  solitons  regain  their  original  amplitudes  and 
speeds  after  a collision;  the  final  effect  of  the  collision 
is  a phase  shift  of  each  soliton.  Two  lumps  regain  their 
original  amplitudes  and  speeds  after  a collision,  and  suffer 
no  phase  shift. 

(iv)  Explicit  formulae  are  available  for  N solitons, 
and  for  N lumps.  The  formulae  for  the  one  soliton  and  one 
lump  solutions  of  (1.3),  with  o = -1  , are: 
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I 


Sol iton 


u 


£.n 


{ 1 + exp (~2rx ' ) } 


x'  ■»  x + 4k  2 1 


4 

(3.1) 


Lump 

u ■=  -12-^r  £n  {(x'  + py')2  + (qy’)2  + 3/q2}  , (3.2) 

u X 

where 

x’  - x + (p2  + q 2 )t 
y'  ■=  y - 2pt 

These  two  solutions  are  drawn  in  Figure  2 for  a particular 
choice  of  the  constants.  (The  soliton  is  a negative  wave 
in  Region  B,  as  shown.  In  Region  A,  solitons  are  positive.) 

These  stability  results  suggest  that  whereas  the  one- 
dimensional KdV  solution  may  play  an  important  role  in 
(1.3)  with  a - +1  (Region  A),  no  such  situation  is  en- 
visaged when  o - -1  (Region  B) . 


I 

I 

I 

I 

I 


3 . 2 The  Nonlinear  Schrodinger  Limit 

Next,  we  consider  the  nonlinear  ' chrodinger  equation  (2.25). 
Observe  that  (2.25)  admits  one-dimensional  solitons  traveling  at 
almost  any  acute  angle  relative  to  the  group  velocity  of  the 
packet.  The  extreme  cases  are  found  by  setting  either  3/3n«  0 
or  t)/3£  ■ 0 . If  3/3n  ■ 0 , the  second  equation  in  (2.25) 
ran  be  integrated  once,  and  the  system  reduces  to 


iAT  + AAU  - v | A | 2A 
- -B/a | A | 2 


(3.3) 
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| 

I 


where  v = x~Xi B/a  , and  the  coefficients  A , a , B , x . 

X:  are  defined  in  (2.26).  (Throughout  this  discussion, 

it  should  be  borne  in  mind  that  the  amplitude  A represents 
the  envelope  of  a train  of  plane  waves.)  Initial  data 
can  be  created  experimentally  by  modulating  (in  time)  the 
stroke  of  an  oscillating  paddle  at  the  end  of  a one-dimensional 
wave  tank.  If  Av  > 0 , as  it  is  in  Regions  A,  B and  E, 
of  Figure  1,  there  are  no  solitons.  The  initial  data 
evolve  into  a field  of  decaying  oscillations  that  we  shall 
refer  to  as  "radiation."  This  radiation  can  be  described 
in  terms  of  a modulated  similarity  solution,  and  it  decays 
as  t~%  [23].  In  Regions  C,  D and  F,  Av  < 0 , and  the 

same  initial  data  now  produce  a finite  set  of  envelope 
solitons  in  addition  to  the  radiation.  (For  appropriate 
initial  data,  multi-soliton  states  are  also  possible  [1]). 

The  one-soliton  solution  of  (3.3)  is 

A = a| 2A/v |^sech{a(^-2bx) }exp{ib^  F iA(a2-b2)i)  . (3.4) 

The  constant  b in  (3.4)  represents  an  0(c)  correction 
to  the  basic  wavenumber,  k ; without  loss  of  generality 
we  take  b = 0 . It  is  evident  from  (3.4)  that  the  am- 
plitude of  the  envelope  soliton  is  of  permanent  form,  and 
represents  a physical  wave  that  decays  only  due  to  viscous 
effects.  Figure  3 shows  the  experimental  measurements  of 
such  a wave,  and  we  have  superposed  on  the  measurements 
the  soliton  solution  with  the  same  peak  amplitude.  [This 
experiment  was  conducted  by  Professor  J.L.  Hammack  while 


L. 
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at  the  University  of  Florida,  and  we  are  grateful  to  him 
for  allowing  us  to  use  his  unpublished  data.]  It  is  clear 
from  this  comparison  that,  at  least  in  some  aspects,  the 
model  represented  by  (3.3)  is  remarkably  accurate.  For 
more  detailed  comparisons,  see  [4]  or  [24]. 

At  the  other  extreme,  if  9/9£  = 0 in  (2.25),  the  system 
reduces  to 

iAx  + uArin  -=  x|A|2A  , (3.5) 

which  is  mathematically  equivalent  to  (3.3)  but  represents 
a much  different  situation  physically.  Here  wave  crests 
move  in  the  x (or  O direction,  but  they  are  modulated 
in  the  q direction.  These  modulations  can  move  only 
in  the  q direction.  To  our  knowledge,  this  configura- 
tion has  not  been  explored  experimentally  in  water  waves, 
although  it  is  common  in  nonlinear  optics,  where  Arr) 
represents  diffraction  of  the  light.  In  optics,  initial 
data  is  produced  experimentally  with  a diffraction  grating, 
and  the  solution  of  (3.5)  provides  a nonlinear  description 
of  Fraunhofer  diffraction  (cf.,  Manakov,  [25]).  Solitons 
exist  where  x < 0 (since  p > 0)  in  Regions  B,  C and  F. 

To  distinguish  them  from  the  soliton  solutions  of  (3.3), 
we  will  refer  to  the  solitons  in  (3.3)  as  "envelope  solitons," 
and  the  solitons  in  (3.5)  as  "waveguides." 

Between  these  two  extremes,  9/9q  - 0 and  9/9£  « 0 , is  a 
one-parameter  family  of  other  one-dimensional  restrictions  of 
(2.25),  corresponding  to  one-dimensional  waves  (of  the  envelope) 
traveling  at  various  angles  relative  to  the  group  velocity  of 
the  carrier  wave.  Each  of  these  one -dimensional  problems  is 
governed  by  an  equation  of  the  form  (1.2),  except  at  one  angle 
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that  corresponds  to  crossing  from  region  B to  F,  and  another 
that  corresponds  to  crossing  from  F to  D. 

Again,  the  question  arises  of  the  physical  relevance 
of  the  one-dimensional  soliton  in  the  two-dimensional 
problem.  For  the  nonlinear  Schrodinger  equation,  (2.25), 
the  answer  seems  to  be  that  except  for  specially  contrived 
one-dimensional  geometries  (like  laboratory  wave  tanks), 
they  are  unlikely  to  persist.  We  show  next  that  every 
one- dimensional  soliton  solution  of  (2.25),  envelope 
soliton  or  waveguide,  is  unstable  with  respect  to  a long- 
wave transverse  perturbation.  Apparently,  this  instability 
has  not  been  observed  in  wave  tanks  only  because  the  tanks 
are  too  narrow  to  admit  the  long-wave  perturbations  re- 
quired. The  instability  was  discovered  first  by  Zakharov 
and  Rubenchik  [12]  for  (2.27).  Our  analysis  is  a generaliza- 
tion of  theirs  to  the  case  of  finite  depth. 

Consider  first  the  envelope  solitons,  which  are  solu- 
tions of  (3.3)  and  can  exist  in  Regions  C,  D and  F in 
Figure  1.  As  remarked  above,  it  is  sufficient  to  demonstrate 
the  instability  of  the  stationary  soliton: 


A = exp(iAa2x)MO 
- -B/aiJ/2(0 

where  MO  is  real  and  satisfies 


(3.6) 


+ Aa2  + vM  - 0 


(3.7) 


Perturbations  about  this  soliton  can  be  put  in  the  form 
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A = exp(iAa2T)[iJ>  + u + iv] 

oo 

<J>  = (B/cx)  f (4*2  + 2ijm)dz  + w , 


(3.8) 


where  u , v , and  w are  real , 

I u | , | v|  <<  <J>  , | w | <<  4>  , 

and 

u,v,w  ~ exp(ipn  ± iftx) 

The  question  of  stability  now  comes  down  to  determining 
whether  ft2  is  positive.  Substituting  (3.8)  into  (2.25), 
linearizing  and  eliminating  v yields 

N 

ft2u  = (Lo  + yp2  ) (L)  + yp2)u  + Xt(L0  + yp2)^w^  , 

OO 

aw^  = P2|^f  J(^u)dz  + wj 


where  Lq  and  Lj  are  the  self-adjoint 
operators  defined  by 


(3.9) 


-X  -11 

+ Xa 

32 

+ Xa 

In  the  short-wave  limit  (p2  -*•«>)  , (3.9)  reduces  to 

ft2u  = y2pl,u  + 0(p2 ) , 

oo 

^ f (<|m)dz  + w ~ 0 


(3.10) 
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Clearly  ft2  is  positive  in  this  limit,  and  short  waves 
are  not  unstable.  Indeed,  if  they  were  unstable,  it  would 
be  difficult  to  observe  envelope  solitons  even  in  narrow 
wave  tanks . 

In  order  to  analyze  the  long  wave  limit  (p2  -*■  0)  , we 
expand  the  unknowns  in  (3.9)  as 

u ~ uo  + p2ui 

w ~ p 2 w ! , (3. 11) 

ft2  ~ p2ft? 

Then  to  leading  order,  (3.9)  becomes 

L L, u = 0 . (3.12) 

o 1 o 

In  order  to  solve  (3.12),  we  define  certain  odd  (-)  and 
even  (+)  functions  of  £ : 

uo  = h ’ Uo  = 4 • vo  = ‘^/2A  • vo  = * • (3-13) 


The  following  relations  can  be  obtained  from  (3.7) 


It  follows  that  u. 


L 

= 0 

o o 

L v 

= u 

o o 

o 

L.  u” 

= 0 

i o 
. + 

+ 

Liuo 

= V 

o 

4* 

uo 

both 

(3.14) 


that  v~  and  v^  satisfy  the  adjoint  equation, 

L , L v = 

1 o o 


0 


(3.15) 
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In  each  case,  there  are  two  other  solutions  that  do  not 
vanish  as  | £ | •+  °°  . We  will  also  need  certain  scalar 
products  of  these  functions.  Using  the  notation 


oo 

(4>  1 , <t>2  ) = 

— 00 

one  computes 


/ 


<P  1 <P  2 d£ 


where 


<V%>  - 
. + +*  1 dl 


<vo’uo>  ■ 

1 dl 

A da2  ’ 

(3.16) 

<v>;>  - 

<Vuo>  = 0 • 

<Vuo>  ■ 

$ 4^2d£  , 

I = J | A | 2 d£  = 4 


At  0(p2)  , (3.9)  reduces  to 

LqL i u j = fi2uo  - u (Lq  + h ! ) Uq  - X i (Wj 

- (3.17) 

“(wi)f^  = ~ / (^uQ)dz 

C 

For  u,  to  decay  as  | £ | -+  » , it  is  necessary  that  the 
nonhomogeneous  terms  in  (3.17)  be  orthogonal  to  the  decaying 
solutions  of  the  homogeneous  adjoint  equation  (3.15). 

Because  the  equations  are  linear,  it  is  sufficient  to  con- 
sider the  odd  and  even  modes  separately.  Thus,  if  uq 
in  (3.17)  is  u+  , we  multiply  (3.17a)  by  v+  , integrate 
over  £ , and  use  integration  by  parts  to  obtain 
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or 


0 * <ai>+<4  ]£■>  - ul 


(Si?)' 


= -2AMa2 
dl 


& 


(3.18) 


For  the  odd  mode,  uo  , we  multiply  (3.17a)  by  v”  and 
use  (3.17b).  The  result  is 


<»?>-  - x[u  / *ldc,  + *'d«] 

(3.19) 

= ^ Aa2|]  + |X/.,] 


The  question  of  stability  of  envelope  solitons  depends 
only  on  the  sign  of  A (the  other  factors  in  (3.18)  and 
(3.19)  are  intrinsically  positive).  Using  ft2  ~ p2(tt2) 
we  summarize  the  result  as  follows: 

(i)  In  Region  D,  where  A < 0 , an  envelope  soliton 
with  amplitude  a is  unstable  with  respect  to  long  dis- 
turbances that  are  antisymmetric  (-)  in  £ . The  growth 

rate  (ift)  of  the  disturbance  with  wavenumber  p i , 
found  from 


02  = -£  p 2 a 2 | A | ( u + | A/ v |)  + 0 (p  “)  (3.20) 

In  the  deep  water  limit,  this  simplifies  to 

02  = -£  p2a2p  | A | + 0(p‘‘)  , (3.21) 

as  found  by  Zakharov  and  Rubenchik  [12].  Thus,  for  an 
inviscid  fluid,  the  effect  of  finite  depth  is  to  enhance 
the  growth  rate  of  the  instability.  Zakharov  and  Rubenchik 
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found  the  OCp1*)  correction  to  (3.21),  and  argued  qualita- 
tively that  the  most  unstable  wave  satisfies 

up2  «=  0 ( | A | a2  ) (3.22a) 

and  that  the  maximum  growth  rate  is  on  the  order  of 

|«|  - 0 ( | A | a2 ) . (3.22b) 

Moreover,  they  noted  that  the  growth  of  a mode  that  is  anti- 
symmetric in  £ and  sinusoidal  in  > tends  to  bend  the 
wavecrest,  producing  a "snake"  effect;  i.e.,  the  crest 
of  the  perturbed  wave  oscillates  back  and  forth  in  the 
(£-n)  plane  about  its  unperturbed  position.  Recent  numerical 
computations  in  [26]  have  made  (3.22)  more  precise. 

(ii)  In  Regions  C and  F,  where  A > 0 , an  envelope 
soliton  with  amplitude  a is  unstable  with  respect  to 
long  symmetric  (+)  disturbances.  The  growth  rate  of  the 
disturbances  with  wavenumber  p is  found  from 

ft2  - -2p2a2 Ap  + O(p-)  . (3.24) 

and  this  result  also  holds  in  the  deep  water  limit.  Again, 
qualitative  considerations  yield  (3.22b).  Growth  of  a 
symmetric  mode  tends  to  modulate  the  wave  amplitude  periodi- 
cally in  n . 

Analysis  of  the  stability  of  waveguides  (in  Regions  B, 

C and  F)  follows  similar  lines,  and  it  is  necessary  only 
to  indicate  the  main  points  of  the  analysis.  A stationary 
waveguide  has  the  form 
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A = exp(iua2T)tKn)  ( 3 . 23 ) 

* « 0 

where  <1/  is  real  and 

-u4'nn  + n a 2 - x'P 3 * 0 . 

Perturbations  take  the  form 


A ■ exp(iua2T)[^  + u + iv]  , 


oo 


(3.26  ) 


The  linearized  equations  for  u and  w , derived  from  (2.25), 
are 


ft2u  * (Lq  + Ap2)(L1  + Ap2)u  + X1(LC  + A p 2 ) vp  w 
wr  r ■ p2 [ 2Blu  + aw]  , 

where 

Lo  “ h7  + ua 2 + x,j;2  * 

- -u  ^rr  + pa2  + 3y>2  • 


■\ 


) (3.27) 


These  equations  are  very  similar  to  those  in  (3.9) 
and  we  simply  state  the  final  result.  Throughout  Regions  B, 
C and  F,  (A,u)  are  positive.  Anywhere  in  these  regions, 
a stationary  waveguide  with  amplitude  a is  unstable  with 
respect  to  long  symmetric  (in  n)  disturbances.  The  growth 
rate  (ift)  of  the  disturbance  with  wavenumber  p (in  O 
is  found  from 
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U2  = -2p2  a2 Ay  + 0(p4)  , (3.28) 

and  qualitative  considerations  given  (3.23a,b)  with  X , p 
interchanged . 

We  conclude  this  section  by  summarizing  our  results 
for  the  nonlinear  Schrodinger  equation,  (2.25).  There  are  many 
one-dimensional  limits,  including  (3.3)  and  (3.5).  These  two  limits 
admit  envelope  solitons  and  waveguides,  respectively,  in 
various  regions  of  Figure  1.  However,  all  possible  solitons 
are  unstable  with  respect  to  some  long-wave  transverse 
perturbation.  This  instability  does  not  appear  in  experi- 
ments in  one-dimensional  wave  tanks,  provided  the  tank 
width  is  small  in  comparison  with  the  soliton  length, 
because  the  unstable  modes  are  excluded  by  the  geometry. 

If  this  constraint  is  removed,  however,  the  instability 
should  occur,  and  neither  kind  of  soliton  is  a stable 
asymptotic  state  that  can  be  achieved  from  initial  data 
in  (2.25). 


4. 


FOCUSING 


In  one-dimensional  problems,  like  (1.2),  the  most  dramatic 
nonlinear  effect  is  that  smooth  initial  data  can  "focus" 
into  a localized  soliton,  or  into  a set  of  solitons,  which 
then  persist  forever.  In  this  section,  we  show  that  focusing 
is  even  more  dramatic  in  two  dimensions  and  that  a solution 
of  (2.25)  that  evolves  from  smooth  initial  data  can  become 
singular  at  a point  in  space  after  a finite  time.  This  is 
known  as  the  "self-focusing  singularity,"  or  simply  as 
"focusing."  In  such  a case  the  water  wave  equations  must  be 
re-examined  in  the  neighborhood  of  the  focus. 

To  our  knowledge,  the  phenomenon  of  focusing  has  not 
yet  been  observed  as  such  in  water  waves,  although  it  has 
been  known  for  some  time  in  nonlinear  optics  (e.g.,  Vlasov, 
Petrishchev  and  Talonov  [27]).  Some  of  the  analysis  dis- 
cussed here  uses  the  ideas  presented  by  Zakharov  and 
Synakh  [28]  who  studied  what  amounts  to  the  two-dimensional 
version  of  (1.2)  (i.e.,  2.27)  in  the  context  of  the  optics 
problem. 

4 . 1 Necessity  of  Focusing 

Our  first  objective  is  to  identify  circumstances  under 
which  the  solution  of  (2.25)  must  focus  in  a finite  time. 
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Consider  any  point  in  Region  F of  Figure  1,  where 
(A , u , (-x) , X j . « » B)  are  all  positive;  i.e.,  consider  capillary- 
type  waves  in  sufficiently  deep  water.  Consider  initial 
data  for  (2.25)  which  are  infinitely  differentiable  and 
which  decay  rapidly  as  (£2  + n 2 ) -*■  <»  ; e.g.,  A(£,n,0) 

might  have  compact  support.  If  a solution  of  (2.25) 
exists  and  vanishes  rapidly  enough  as  (£2  + n2)  -*•  °°  , 
then  the  following  integrals  are  constants  of  the  motion: 


Ij_  = ff  | A | 2d£dn 


h '//  <A  TC  - A*  W)iKin 


Each  bracket,  {}  , in  1^  is  positive  definite,  and 
the  second  bracket  vanishes  in  the  linear  limit  of  (2.25). 
Clearly  1^  <0  is  possible  (e.g.,  if  the  initial 
data  has  sufficiently  large  amplitude) . 

It  also  follows  from  (2.25)  that 


If  {x  + t}  |A|2dcdn  " 8I4 


(4.2) 
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As  noted  in  Section  2,  one  may  interpret  1^  as  the  mass 
of  the  wave  (to  leading  order  in  e)  . Then  the  integral 
in  (A. 2)  may  be  interpreted  as  the  moment  of  inertia,  and 
(A. 2)  is  an  example  of  the  virial  theorem  (e.g., 

[29],  P-  581).  (A. 2)  is  easily  integrated,  and  we  see  that 

if  1^  < 0 , then  the  moment  of  inertia  vanishes  at  a finite 
time.  Clearly,  no  global  solution  exists  after  this  time, 
because  the  (positive  definite)  moment  of  inertia  would 
become  negative!  Since  the  mass  of  the  wave  is  conserved, 

(A. 2)  suggests  that  prior  to  this  time  the  radius  of  gyration 
is  vanishing  as  the  mass  accumulates  at  a single  point.  The 
rapid  development  of  this  singularity  is  what  we  mean  by 
focusing . 

Before  examining  the  nature  of  the  singularity  that 
develops,  let  us  consider  the  implications  of  this  argu- 
ment outside  of  Region  F.  In  Regions  B and  C,  where 
a < 0 global  integrals  involving  <1  are  generally  un- 
bounded (cf . , Section  2)  and  no  global  information  about 
the  solution  is  available  by  this  approach.  Whether  focusing 
exists  in  these  regions  is  open.  In  Region  E there  is  no 
focusing  in  the  deep  water  limit,  since  the  parameters  are 
such  that  1^  > 0 . In  arbitrary  depth  the  question  of 
focusing  is  still  open. 

In  Regions  A and  D,  the  integral  in  (A.  2)  is  not  of 
definite  sign,  and  provides  no  contradiction.  Both  be- 
cause of  the  breakdown  of  this  argument  and  because  the 
type  of  instability  of  solitons  is  different  than  in  Region  F, 

I 
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we  expect  that  if  singularities  develop  in  these  re- 
gions, they  will  be  qualitatively  different  than  those 
of  the  self-focusing  type. 

4 . 2 Nature  of  the  Singularity 

Next,  we  examine  the  possible  behavior  of  the  singularity 
that  develops  at  the  focus.  Zakharov  and  Synakh  [28] 
studied  the  radially  symmetric  case  of  (2.27).  They  in- 
vestigated this  equation  both  by  numerical  compuations 
and  an  approximate  analytic  procedure.  From  these  they 
concluded  that  as  t •+  tq  (tq  being  the  time  of  focus) 
the  wave  amplitude  grows  as  (t  -t)“P  , p *=  2/3  . 

In  this  section  we  show  that  there  are  a number  of 
quasi-self-similar  solutions  to  the  generalized  nonlinear 
Schrodinger  equation,  (1.4),  including  one  with  p * 2/3  , 
but  we  have  found  no  convincing  argument  that  this  local 
behavior  is  necessarily  of  the  p “ 2/3  type. 

For  convenience,  we  consider  the  scaled  form  of  (2.25), 
namely  (1.4).  In  Region  F of  Figure  1,  where  focusing 
can  occur,  0j  ■=  +1  , o2  *=  -1  . 

Let  A ■ B exp(i'i')  in  (1.4),  with  B,Y  real  and 

find  : 

(%B2),+  (*  B2)  + (y  B2)  - 0 , (4.3a) 

L x x y y 

B + B + B - B(y2  + y2)  - -B3  + # B 
t xx  yy  VTx  y7  *x 


(4.3b) 
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a*!1  + 4>  *=  -b(B2)  (A. 3c) 

xx  yy  v ' x v ' 

We  seek  quasi-self-similar  solutions  of  (4.3)  in  the  neigh- 
borhood of  the  point  of  focus  in  the  form 


- ^ R(x,y)  + RQ(x,y,t) 

(4.4a) 

- j Q(x.y)  + Qo(x,y,t) 

(4.4b) 

where  x = x/f  , y = y/f  , f(t)  = (tQ-t)^  , so  that  f 0 
as  t -*■  tQ  . This  expansion  is  asymptotic  near  the  focus  provided 
RQ  <<  R . Qq  <<  Q in  this  region.  Zakharov  and  Synakh  [23] 
also  assumed 


R0  - 0(fR)  , (4.5) 

but  this  assumption  seems  to  be  unnecessary.  In  any  case,  the 
dominant  terms  in  (4.3a)  as  t t are 

(^R2  - £ R2f  ')x  + (4>yR2  - | R2f  ')y  - 0 (4.6) 

A special  solution  of  (4.5)  is 

y = £ f'  + 
x ^ R 

(4.7) 

Ty  * \ 

Taking  Gi  - G2  =0  (for  which  some  motivation  is  provided  below) 
' yields 

V m ~7~  (x2  + y2)  + g(t)  , (4.8) 


and  with  this  we  have  from  (4.3)  - (4.4),  as  f -►  0 
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R--  + R-- 
xx  yy 

+ R3  - 

RQ-  - g'( t ) f 2 R - 

t 

- (x2  + y2 )R  ~ 0 (4.9a) 

aQ;x  + 

Q-  - + b (R2  ) - ~ 0 
xyy  x 

• 

(4.9b) 

There  are 

various  possibilities; 

e-g. 

$ 

(a)  f"f3 

<<  1 

. g'=  </f2  ; 

(b)  f"f3 

= 0(1) 

. g’-  </f2  ; 

others  are  obtained  similarly.  Case  (a)  implies  that 
to  leading  order  (4.9)  reduces  to 

R- - + R- - + R3  - RQ-  - tcR  « 0 , (4.10a) 

xx  yy  Xx  ’ v 

aqJ5  + Qyj,  - -b(R>)j  . (4.10b) 

if  one  also  assumes  (4.5),  then  (a)  becomes  f3f"  * 0(f)  , from 
which  it  follows  that  p - 2/3  . However,  the  spatial  structure 
defined  by  (4.10)  does  not  depend  on  (4.5),  or  on  p « 2/3  . 

In  the  deep  water  case  with  radial  symmetry  [23],  b = Q ■=  0 , 

?2  = x2  + y2  , and  (4.10a)  reduces  to 

Rii  + = Ri  + R!  - *R  - 0 • <41l> 

Chiao,  Garmire  and  Townes  [30]  first  studied  (4.11)  as  a 
model  of  cylindrical  optical  beams,  and  showed  that  its 
bounded  solutions  decay  exponentially  for  large  r . The 
equation  also  arises  as  an  exact  reduction  of  (1.4)  if 
we  take 


-40- 


I 

| 

I 

I 

I 

I 

I 

I 

I 


Q - b - 0 

r “ A/x2  + y2  , (4.12) 

A - AR(r)exp(itcA2t)  , 

The  fact  that  (4.11)  is  exact  has  important  consequences, 
which  we  discuss  in  Section  5. 

In  case  (b) , p * 1/2  and  the  solution  is  exactly 
self-similar.  Here  (4.8)  - (4.9)  yield 

R—  — + R- - + R3  - RQ-  - kR  + X (x2  + y2)R  ■=  0 (4.13a) 

xx  yy  16 

aQji  + Qyy  - -b(R!)j  (4.13b) 

y = --g-(x2  + y2)  - Kin  (tQ  - t)  + 4'q  . (4.13c) 

In  the  radially  symmetric  case,  (4.13a)  becomes 

R--  + 1 r-  + - k)r  + R3  = o , (4.14) 

and  for  large  r all  bounded  solutions  decay  as  (r)-1  . 

We  also  note  that  a somewhat  more  general  equation  than 
(4.14),  obtained  by  retaining  G in  (4.7),  can  be  found 
in  the  symmetric  case: 

R?f  + FRf  + <!f-  k>r  + r’  - rw-  0 <A15a) 

r 

, . -|1-  Kt„(to  . t)  +0/^^+  T„  (4.15b) 

However,  one  can  show  from  (4.15a)  that  R has  a finite  value 
at  the  origin  only  if  C « 0 . This  result  provides  some 
justification  for  neglecting  G in  (4.7). 
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Using  any  of  these  similarity  solutions,  going  back  to 
the  full  water  wave  equations  and  rescaling,  we  find  that  the 
focusing  instability  produces  a finite  (i.e.,  0(1))  region 

of  space  in  which  the  wave  amplitudes  are  potentially  large 
enough  to  break  (0(1))  . 

Next,  we  present  an  argument  which  suggests  that  the 
self-focusing  singularity  cannot  be  of  the  p ■ 1/2  type, 
as  we  have  described  it  here.  For  convenience,  we  consider 
the  case  of  radial  symmetry.  The  expansion  in  (4.4)  is  valid 
in  a region  near  the  focus  where  the  first  terms  are  dominant. 
If  we  assume  that  this  "inner  solution"  matches  to  an  outer 
solution  that  is  0(1)  , then  the  inner  expansion  breaks  down 
where 


£ R (?)  “ 0(1)  (4.16) 

But  for  large  ? , R decays  as  (r)_1  , so  that  ^ R 
decays  as  r-1  ; i.e.,  there  is  no  time-dependence.  It 
follows  that  (4.16)  defines  a boundary  for  the  focal  region, 
denoted  by  r - 0(L)  , where  L is  time-independent.  The 
mass  within  this  region  is  proportional  to 

L L/f 

M - f r R2 (r)  dr  - f rR2(r)dr  , (417) 

o 0 

and  {because  R - r_l)  this  grows  logarithmically  as  t -*■  t . 
But  the  total  mass  is  finite,  and  this  is  a contradiction. 

In  case  (a)  R(r)  decays  exponentially  and  no  such 
contradiction  appears.  Moreover,  if  the  nonlinearity  in 
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(1.4)  were  slightly  stronger,  no  contradiction  appears  in 
the  purely  self-similar  case.  To  be  precise,  if  the  non- 
linear term  in  (1.4)  were  replaced  with  | A | 2 aA  , a > 1 , 

then  the  radially  symmetric  similarity  solution  becomes 

_ 1 

A = (tQ-t)  , i « r/(to-t)%  . (4.18) 

- -1/a 

In  this  case  the  radial  decay  is  B(r)  - (r)  , and 

again  L is  finite.  However,  in  this  situation,  when 
a > 1 the  mass  remains  finite  as  t -*•  t , and  the  pure 
similarity  solution  is  a likely  candidate  for  describing 
the  dynamics  of  the  focus  regime. 

Finally,  we  remark  that  a natural  generalization  of 
(1.4)  is: 

iA  + o i A + A + o2A  „ - o3|A|2aA  + * A , (4.19a) 

t 1 xx  yy  1 zz  i 1 1 x ’ 

ai4xx  + *yy  + aj*zz  " -b(lAl2a>x  . (4.19b) 

where  the  ck  ■ ±1  (i  - 1,2,3)  and  a\  , 82  , b are 
constant.  The  spherically  symmetric  limit  is  obtained  by 
taking  b “$“0  , Oj  - o2  ■ +1  . Since  the  spherically 
symmetric  equation  has  wide  applicability,  and  (1.4)  is 
itself  physically  relevant,  we  expect  that  (4.23)  will 
also  arise  in  physical  problems. 
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5.  OTHER  SOLUTIONS  OF  THE  NONLINEAR  SCHRODINGER  EQUATION 

The  purpose  of  this  section  is  to  identify  other  fea- 
tures of  the  solution  of  (2.25)  that  may  play  a role  in 
its  asymptotic  (t  -*■  °°)  solution. 

5 . 1 Complete  Integrability 

Perhaps  the  fundamental  question  to  answer  about  (2.25) 
is  whether  it  is  completely  integrable;  i.e.,  whether  it 
can  be  solved  exactly  by  relating  it  to  an  appropriate 
linear  scattering  problem.  The  question  is  natural  in 
light  of  the  fact  that  the  one -dimensional  problem  can  be 
solved  in  this  way. 

Consider  first  the  long-wave  limit  of  (2.25),  subject 
to  the  constraint  in  (2.16).  Here,  (2.25)  becomes  (after 
rescaling  of  variables) 

iAT  - °iA^  + A^n  = o,  | A | 2A  + A4>£ 

(5  • 

0!^  + »nn  -=  -2  ( I A I 2 )^  , Oj  « sign(^-T)  . 

This  system  is  of  I.S.T.  type  [11].  Special  N solitons 
solutions  can  be  constructed  either  by  a direct  (Hirota 
type)  method  or  via  the  Zakharov-Shabat  approach  [14,31]. 
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The  situation  seems  to  be  much  different  in  the  deep 
water  limit.  Here  we  have  already  seen  that  (4.11)  is  an 
exact  reduction  of  (2.25)  to  an  ordinary  differential 
equation;  i.e.,  every  solution  of  (4.11)  provides  an  exact 
solution  of  (2.25)  in  this  limit.  Let  us  consider  those 
partial  differential  equations  (PDE)  which  have  been  solved  . 
exactly  by  1ST  methods.  We  have  found  that  every  reduc- 
tion of  one  of  these  PDE's  to  an  ordinary  differential 
equation  (ODE)  results  (perhaps  after  a transformation 
of  dependent  variables)  in  an  ODE  without  moveable  critical 
points  [32 , 33  ] . 

We  expect  that  if  (2.25)  can  be  solved  by  some  1ST, 
then  (4.11)  should  have  no  moveable  critical  points.  But 
Ince  [34],  esp.  p.  344)  provides  a complete  list  of  all 
such  second-order  equations;  (4.11)  is  not  on  this  list 
and  cannot  be  transformed  to  an  equation  on  this  list. 
Therefore,  the  solution  of  (4.11)  has  moveable 

critical  points.  Moreover,  one  can  show  that  (4.11) 
has  logarithmic  singularities  in  addition  to  poles.  On 
this  basis,  we  conjecture  that  (2.25)  cannot  be  solved 
exactly  by  1ST  in  the  deep-water  limit. 

Although  (2.25)  can  be  solved  by  1ST  in  the  shallow- 
water  limit  (i.e.,  lower-left  corner  of  Figure  1),  it 
apparently  cannot  be  solved  in  this  way  in  the  deep-water 
limit.  Wherever  1ST  methods  fail,  one  is  forced  to  piece 
together  special  solutions  of  the  problem  to  describe  the 
general  solution. 
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5 . 2 Decaying  Oscillations  1 

The  special  solutions  discussed  so  far  in  this  paper  j 

have  been  localized:  either  solitons  (or  soliton- like) 

or  self-focusing  singular  solutions.  However,  in  the  one- 

dimensional  limit  of  (2.25),  solitons  make  up  only  part 

of  the  asymptotic  solution  of  the  initial  value  problem. 

That  part  of  the  solution  associated  with  the  continuous 

spectrum  spreads  over  large  regions  of  space,  while  it 
-k 

decays  as  t . In  particular,  an  exact  solution  of  • 

(1.2)  is 

A = t~^  A exp{i(^  + oA2  In  t + 0)  } , (5.2) 


where  A and  0 are  real  constants;  the  solution  of 
(1.2)  associated  with  the  continuous  spectrum  tends  to 
a slowly-varying  modulation  of  this,  where  A and  0 
depend  on  (x/t)  [22,35]. 

In  the  two-dimensional  problem,  (1.4),  there  is 
an  analogous  exact  solution: 


A - t_1  A exp{i  (qiX4t+-y  - + a2  A2/t  + B(t)  + 0)  } 
4-  - -B’  ( t ) x + C(t)y  + D(t)  . 


(5.3) 


Similar  solutions  in  the  deep  water  limit  of  (2.25)  were 
found  by  Talanov  [36].  On  the  basis  of  the  one-dimensional 
theory,  we  anticipate  that  the  part  of  the  solution  of  (1.4) 
that  decays  in  time  can  be  described  in  terms  of  a slowly- 
varying  modulation  of  this  exact  solution. 
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Moreover,  this  behavior  would  be  consistent  with  the 
results  of  Lin  and  Strauss  [37]  who  studied  the  three- 
dimensional  problem 

iut  - Au  + |u|2u  =•  0 , (5.4) 

where  A is  the  Laplacian  in  three  dimensions.  They  found 

-3/2 

that  the  solution  exists  for  all  time  and  decays  as  t . 

The  appropriate  similarity  solution  here  is 

u = t"3/2  A exp { i (-*I+Xl±£l)  + A2/2t2  + $)}  . (5.5) 

Without  solitons  or  focusing,  the  decay  rate  of  the  solution 
of  the  nonlinear  Schrodinger  equation  seems  to  be 

u = 0(t‘n/2)  (5.6) 

where  n is  the  number  of  spatial  dimensions.  This  decay 
rate  is  the  same  as  in  the  linearized  problem. 
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LIST  OF  FIGURE  CAPTIONS 

Figure  1.  Map  of  parameter  space,  showing  where  the  coef- 
ficients in  (2.25)  change  sign.  The  dynamics 
of  wave  evolution  is  different  in  each  region. 

Figure  2a.  KdV  soliton,  as  seen  in  two  space  dimensions 

at  a fixed  time;  k2  = 1/12  in  (3.1),  with  o = -1 

Figure  2b.  Lump  solution  of  (3.2)  as  seen  in  two  dimensions 
of  a fixed  time;  p = 0 , q2  = 1/3  , a - -1  . 

Figure  3.  Measured  surface  displacement,  showing  evolution 
of  envelope  soliton  at  two  downstream  locations; 
h = 1 m , kh=4.0,  w = 1 cps  , 

T = 1.0  * 10-“  ; 

, measured  history  of  surface  displacement; 

, theoretical  envelope  shape 

<4  = <a  sech(z) 

z = (V/8A)*(C  t-x) 

0)  N g 

(3a)  6 m downstream  of  wavemaker,  xa  = 0.132  . 
(3b)  30  m downstream  of  wavemaker,  xa  = 0.116  . 

Figure  4.  Stationary  waveguide,  as  seen  in  time  at  a fixed 
location.  In  (3.5),  y = 2 , X ~ and 

Re(A)  = 2 sech  2p  cos  8x  is  plotted.  The  dis- 
placement of  the  free  surface,  x£  , is  similar. 
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Fig.  2b. 
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